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Abstract 

Bubonic plague (a fatal, flea-transmitted disease) remains an international public health concern. Although our 
understanding of the pathogenesis of bubonic plague has improved significantly over the last few decades, researchers 
have still not been able to define the complete set of Y. pestis genes needed for disease or to characterize the mechanisms 
that enable infection. Here, we generated a library of Y. pestis mutants, each lacking one or more of the genes previously 
identified as being up-regulated in vivo. We then screened the library for attenuated virulence in rodent models of bubonic 
plague. Importantly, we tested mutants both individually and using a novel, "per-pool" screening method that we have 
developed. Our data showed that in addition to genes involved in physiological adaption and resistance to the stress 
generated by the host, several previously uncharacterized genes are required for virulence. One of these genes {ympth66c, 
which encodes a putative helicase) has been acquired by horizontal gene transfer. Deletion of ympti .66c reduced Y. pestis' 
ability to spread to the lymph nodes draining the dermal inoculation site - probably because loss of this gene decreased the 
bacteria's ability to survive inside macrophages. Our results suggest that (i) intracellular survival during the early stage of 
infection is important for plague and (ii) horizontal gene transfer was crucial in the acquisition of this ability. 
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Introduction 

Plague (caused by the bacterium Yersinia pestis) is one of tlie 
diseases that has most greatly affected the course of history [1,2] 
and is still an international public health concern. The WHO has 
even categorized this disease as re-emergent in some parts of the 
world. This problem is heightened by the emergence of multidrug- 
resistant strains and the potential use of plague as a bioweapon 
[3,4] . Plague is commonly acquired through the bite of an infected 
flea, when Y. pestis is regurgitated into the dermis during the flea 
blood meal [5]. It is thought that Y. pestis is then taken up by 
macrophages, within which it replicates and produces antiphago- 
cytic factors [6] . Next, bacteria released from dying macrophages 
spread into the draining lymph node, where they actively replicate, 
circumvent phagocytosis by polymorphonuclear neutrophils 
(PMNs) and confront a number of toxic effectors released by 
degranulating or lysing PMNs [7,8]. These events lead to a 
hemorrhagic, edematous, swollen and painful lymph node (the 
bubo), which is the defining clinical feature in the bubonic form of 
plague [9,10]. Overwhelmed, the lymph node allows the bacteria 
to disseminate into the blood stream; this produces a fulminant, 
systemic infection and fatal septicemia. In a very low proportion of 
infections, septicemia causes pneumonia, which in turn enables Y. 



pestis to be transmitted from person to person via contaminated 
aerosols. However, we lack important data on the above- 
mentioned mechanisms. For instance, genes dedicated to the 
intracellular step have not been identified. Furthermore, recent 
work suggests that transit through the flea vector enhances Y. pestis' 
ability to resist phagocytosis by the macrophage [11]. Very few 
pathways for resisting innate immune effectors have been 
elucidated [12]. Lastly, the metabolic pathways that enable Y. 
pestis to replicate during infection are incompletely understood. 

Regurgitation of Y. pestis into the skin does not always cause 
bubonic plague and bubonic plague does not always progress 
consistently to septicemic plague. This indicates that the skin and 
the lymph node are the two decisive arenas in which Y. pestis must 
deploy its arsenal to produce a successful infection. Better 
knowledge of this arsenal will undoubtedly improve our under- 
standing of the pathogenesis of plague and should help us to 
develop better countermeasures against this disease. To this end, 
we previously determined the transcriptome of Y. pestis in the most 
relevant context of infection: the rat bubo [7]. Although the study 
results were valuable, they were not predictive of the requirements 
for virulence and merely provided clues about putative biological 
roles of upregulated genes in vivo that will have to be confirmed 
using traditional experimental approaches. We therefore pursued 
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Author Summary 

In order to understand and combat infectious diseases, it is 
essential to characterize the full set of genes required by 
pathogenic bacteria to overcome the many immunological 
and physiological challenges encountered during infec- 
tion. Here, we used a genome-scale approach to identify 
genes required by the bacterium Yersinia pestis in the 
production of bubonic plague (a fatal, flea-borne zoono- 
sis). Our results suggest that when colonizing the 
mammalian host, the bacterium (i) relies on carbohydrates 
as its carbon source, (ii) shifts to anaerobic respiration or 
fermentation and (iii) experiences and resists several (but 
not all) types of stress generated by the host's innate 
immune system. Strikingly, only a small set of genes 
(including horizontally acquired and uncharacterized 
sequences) are required for these infectious processes. 
Further investigations of the ypmtl,66c gene provided 
evidence to suggest that accretion of genetic material via 
horizontal transfer has played a key role in Yersinia pestis' 
ability to successfully initiate infection after the dermal 
fleabite. Lastly, we believe that (i) application of our 
approach to other pathogens and (ii) additional studies of 
selected Yersinia pestis genes important for plague 
pathogenesis (some of which are shared with other 
pathogens) will provide a better understanding of bacte- 
rial pathogenesis in general. 

our efforts by devising a means of screening a mutant T. pestis 
library in which genes previously identified as being up-regulated 
in vivo had been deleted. This method enabled us to draw up a 
hierarchical table ranking genes in order of their importance in 
disease production. The results of further investigations suggest 
that horizontal gene acquisition had a key role in establishing T. 
pestis' ability to successfully initiate infection after the dermal 
fleabite. 

Methods 

Ethics statement 

Animals were housed at the Institut Pasteur de LUle animal 
facility (accredited by the French authorities for the performance 
of experiments on live rodents (#B59-350009) and which complies 
with French and European regulations on the care and protection 
of laboratory animals (EC Directive 86/ 609 and the French Act 
#2001-486, issued on June 6, 2001)). At the time some 
experiments were performed, regulations did not require ethical 
approval. However, all animal experiments had been authorized 
by the French veterinary authorities (#59-350218) and were 
performed in compliance with the terms of NIH Animal Welfare 
Assurance #A5476-01 (issued on July 2, 2007). Since 2013, 
animal experiments were approved by the Nord-Pas-de-Calais 
Ethics board committee (#CEEA 23012). Human serum samples 
were obtained after the provision of written, informed consent by 
the donors. The use of human samples has been declared and 
approved by the MESR and CPP boards (#DC-2013-1765). 

The strains, the mutant library and the plasmids 

The bacterial strains and plasmids used in this study are 
presented in Table S 1 and Text S 1 . Mutant T. pestis strains were 
generated using allelic exchange methods (based on lambda Red 
recombinase or the pCVD442 suicide vector [13,14]) and with the 
primer sets and plasmids described in Tables SI and S2. The 
coding sequence of interest (apart from the first and last 50 base 
pairs) was replaced by a kanamycin, zeocin or trimethoprim 



resistance gene that had been amplified from the suicide vectors 
listed in Table SI. Each mutant's genotype was confirmed by PGR 
assays using primer sets that were complementary to the target 
flanking sequence not involved in the allelic exchange (Table S2). 
Mutant strains were complemented with the recombinant pCRII 
plasmid (Life technologies. Saint Aubin, France) containing a wild- 
type (WT) copy of the gene of interest (including the latter's 
putative promoter sequence). 

Determination of bacterial virulence in rodent models of 
bubonic plague 

Virulence attenuation was evaluated in either pools of mutants 
or individual mutants. Wild-type and mutant strains were cultured 
in Luria broth (LB) at 21°C, quantified (via the optical density at 
600 nm), diluted to the desired cell density in sterile phosphate- 
buffered saline (PBS) and mixed with equal colony-forming unit 
(CFU) counts of other mutants if required. Fifty microliters of the 
bacterial suspension were intradermaUy inoculated (as previously 
described [9]) into 8-week-oId female Brown-Norway rats (Janvier, 
France) or OF-1 female mice (Charles River, France). For 
individual screening, rats and mice were inoculated with 20 and 
10 CFUs respectively. For per-pool screening, rats were injected 
with a total of 100 CFUs (20 CFUs of each of five mutants). 
Groups of eight and ten animals were respectively used for 
individual and per-pool screening experiments. For per-pool 
screening, animals were euthanized when signs of terminal plague 
appeared. The inguinal lymph nodes proximal to the inoculation 
site and the spleens were immediately collected after euthanasia 
and then triturated. Serial dilutions of the triturated organs were 
plated on blood agar containing the antibiotics of interest. The 
CFUs were counted after a 48-hour incubation at 28°C. The 
virulence attenuation of T. pestis was measured by calculating the 
average relative competitive index (ARCI), as follows: 

ARCIx = '- where x is the CFU count of the tested 

n 

mutant, Xo is the count for the most abundant mutant (i.e. the 
strain that acts like a WT strain in the animal of interest), / is the 
animal of interest, n is the number of animals used in the 
experiment, and x<Xo. Groups of five 8-week-old female Brown- 
Norway rats (inoculated as above) were used to study the time 
course of lymph node, blood and spleen colonization by bacteria. 
Lastly, in order to study the role of PMNs cells during infection, 
groups of 8-week-old OF- 1 mice were inoculated intraperitoneaUy 
with 100 |Xg of anti-GRlH- (RB6-8C5) antibody (as described 
above) on the day before the intradermal challenge [8]. 

Growth in human serum 

T. pestis grown overnight in LB at 21°C were centrifuged, 
washed and adjusted to 2x10 bacteria/mL in PBS. Twenty 
microliters of the bacterial suspension were added to 180 |J.l of 
fresh human serum in the wells of a 96-weU plate. After different 
incubation times at 37°C, serial dilutions of the serum were plated 
on blood agar and the CFUs were subsequently quantified. 

Bacterial survival in macrophages 

Macrophages were cultured in high-glucose (4.5 g/L) Dulbec- 
co's Modified Eagle Medium (DMEM) supplemented with 5% 
heat-inactivated fetal bovine serum. All incubations were per- 
formed at 37°C in a 5% COj atmosphere. In the intracellular 
survival assay, cells were infected as previously described [15]. The 
day before interaction with T. pestis, 1.4x10'^ RAW macrophages 
in the incubation medium were allowed to seed in the wells of a 
24-well plate. The macrophages were washed three times in 
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incubation medium prior to addition of T. pestis that had been 
grown overnight at 21°G in LB, washed and suspended in 
DMEM. The multiplicity of infection was 10. Immediately after 
mixture of the bacterial suspension and the macrophages, the 
plates were centrifuged (at 18 g for 5 minutes) and then incubated. 
After 30 minutes, the cells were washed three times. Fresh 
medium supplemented with 8 ng/ mL gentamicin was then added. 
One hour later, the wells were washed twice and fresh medium 
containing 2 (ig/mL of gentamicin was added. At various time 
points after contact, macrophages were lysed by incubation in 
500 )xl of cold water for 10 minutes on ice. The lysate was plated 
on blood agar and the number of CFUs was counted. 

Statistical analysis 

A mutant was considered to be virulence-affected when the 
survival curves for animals infected with the mutant and the wild- 
type strain were found to be significantly different (p<0.05) in a 
Gehan-Breslow-Wilcoxon test. Two-way analysis of variance was 
also used to compare data obtained from in vitro experiments; a p< 
0.05 was considered significant. 

Results 

More than 15% of the Y. pestis loci up-regulated in vivo 
are required for virulence 

We previously reported that 526 T. pestis genes (other than 
previously characterized virulence factors) were up-regulated more 
than two-fold in the rat bubo, relative to T. pestis cultures grown at 
21 or 37°C [7]. Seventy-nine of these are highly likely to encode 
factors important for plague pathogenesis because they are (i) up- 
regulated in both the rat bubo and the mouse lung or (ii) thought to 
be involved in the bacterium's defensive response against nitric 
oxide, oxidative stress and iron stress generated by the innnuiic 
system during bubonic plague in rat [7,16]. We therefore generated 
a library of T. pestis mutants in which one or more of these 79 genes 
had been deleted and then tested the individual mutant's virulence 
in mouse or rat models of bubonic plague. In order to reduce the 
number of mutants and the number of animals to be sacrificed, we 
deleted blocks of neighboring genes (whenever possible and 
regardless of their genetic organization and relationship). Hence, 
a total of 33 mutants (rather than 79) were produced. The mutant 
lacking the functional pyruvate dehydrogenase AceEF [17] was the 
only one to show a serious growth defect on LB agar. As expected, 
this mutant did not kiU any mice 1 5 days after inoculation with a 
dose (-20 CFUs) that usuaUy kills all the animals (Table S3) - 
indicating that the mutant was virulence-attenuated. Screening of 
the library for virulence revealed five attenuated mutants (Table S3). 
The strains lacking the homoserine O-succinyltransferase gene metA, 
the ribonucleoside reductase operon nrdHIEF or the uncharacter- 
ized gene ypo0426 required significantly more time to cause fatal 
plague (median survival a 5 days vs 4 days for the mutants and the 
wild-type strain respectively). The lSypo0988, AptsG and AznuABC 
ypo2062 mutants (respectively lacking the genes for an uncharacter- 
ized protein, the glucose transport PtsG and both the Zinc (Zn) 
transporter ZnuABC and the putative zinc murein DD-endopep- 
tidase YebA [18]) caused plague but at a lower incidence (survival 
rate >60"/(} vs £ 1 2.5 for the mutants and the wild-t}'pe respectively). 
Further virulence assays showed that deletion oiyebA (but not znti) 
was responsible for the loss of virulence in the AznuABC ypo2062 
mutant (Figure 1). 

Next, we sought to establish which of the remaining 447 up- 
regulated genes were required for plague progression in rats. To 
this end, we used the approach described above to build a library 
of 137 mutants and designed a per-pool screening method for 



application in the most relevant model of infection. Hence, rats 
were intradermaUy inoculated with about the dose of T. pestis (the 
input pool) that is regurgitated by the rat flea Xenopsylla cheopis (i.e 
100 CFUs) [19]. In order to reduce the selection of false positives 
and false negatives, each mutant in the pool was given at the 
lowest dose that consistendy killed all the rats in experiments with 
the WT (20 CFUs). In other words, each pool comprised five 
mutants. A wild-type strain was not included in each pool because 

(i) we postulated that at least one mutant in the pool would ha\c a 
wild-type like phenot)'pe (i.e. colonizing and killing the animals as 
efficiently as the wild type) and (ii) \v(; expected a bottleneck 
restricting bacterial dissemination from the skin to the other tissues 
[20]. To overcome bias due to bottlenecking, we were therefore 
obhged to evaluate each pool in several animals (see below). To 
identify mutants within a pool, we provided each one with a 
specific antibiotic resistance marker (kanamycin, trimethoprim or 
zeocin, denoted respectively as K, T and Z). The mutants were 
generated using a virulent strain bearing (or not) a zeocin 
resistance gene at the only neutral att Tn7 chromosomal site 
(Figure SI). Hence, five groups of tagged mutants were produced 
(K, T, Z, ZK, and ZT). Each mutant was counted in the lymph 
nodes and the spleen (the output pool), in order to (i) identify genes 
that might be necessary during early or late stages of infection and 

(ii) reduce the likelihood of missing false negatives (i.e. attenuated 
mutants that might artificially appear to be virulent as a result of 
transcomplementation of the mutant's defects by other strains in 
the pool or the production of an abnormal environment that 
facilitates the growth of an otherwise virulence-attenuated mutant. 
We expected rescued mutants to be more outcompeted in the 
spleen than in the lymph node. Hence, we did not expect a 
virulence-attenuated mutant to be fuUy rescued in the final 
"bacterial read" in the spleen - even when it may have been 
transcomplemented in the pool. Lastly, the bacterial loads for each 
mutant were determined in animals displaying terminal clinical 
signs of plague; at this stage of the disease, (i) the bacterial load 
recovered from an organ (>10' CFU) is high enough to ensure 
reproducible recovery of each mutant, and (ii) all the mutants 
would have experienced similar in vivo growth conditions because 
the organs from terminally iU animals described in the literature 
(and certainly from our experience) were all histologii:ally similar 
[9,21]. Using this strategy, a total of 24 pools comprising five 
mutants and one pool comprising four mutants (i.e. 124 mutants in 
total) was tested. Due to technical issues, the remaining 13 mutants 
had to be individually evaluated for virulence. One of the 13 had a 
longer incubation period in animals and was found to lack the 
prophage j/)o7('fii'^/»o 1091 locus (Table S3). 

Terminal signs of plague occurred 2 to 5 days after inoculation 
of the pools. The mean time to sacrifice was 2.6±0.4 days (Figure 
S2). At the time of sacrifice, the lymph nodes and the spleens from 
all rats were heavily colonized, with a median bacterial count of > 
lO** CFU (Figure S3). These data suggested that (-ach pool 
behaved collectively like a WT strain and that at least one mutant 
within each pool was fuUy virulent because i) the time course of 
infection of animals infected with either ~20 or ~100 wild-type 
CFU did not differ significantly (data not shown) and ii) organs 
from animals showing terminal signs of plague reportedly contain 
more than lO' CFU [9,22,23]. Analysis of the CFU counts in the 
output pools revealed that there was a bottleneck that affected 
bacterial dissemination from the skin to the other tissues. The 
proportions of clones recovered from the lymph nodes and the 
spleens did not appear to follow a particular trend; a clone that 
completely dominated colonization in one animal could be totally 
outcompeted in another (Figure S4). This variable "dominant 
effect" hampers the identification of genes needed for microbial 
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Figure 1. Incidence of plague in mice infected with K pestis lacking the zinc transporter Znu and/or the putative zinc murein DD- 
endopeptidase YebA (YPO2062). (A), mice were injected intradermally with —10 Y. pestis WT (winite circles) and AznuABC ypo2062 (black squares) 
or AznuC (black diamonds). (B), mice were inoculated with as in (A) with WT (white circles) or Aypo2062 (grey circles) mutants. Only the loss of 
YPO2062 significantly decreased the virulence of Y. pestis (p<0.05, relative to WT). Groups of 8 mice were used. Data obtained with the AznuABC 
ypo2062 mutant are representative of two experiments. Data obtained with the AznuC and Aypo2062 mutants come from one experiment each. 
doi:1 0.1 371 /journal.ppat.1 004029.g001 



pathogenesis [24,25]. However, we considered that the "dominant 
effect" observed in individual animals would be balanced by an 
analysis of the group of animals (notably by averaging the data 
obtained from the different rephcates). Hence, we determined the 
relative virulence of each mutant for a specific organ by calculating 
an average relative competitive index (ARCI) for each mutant by 
using the equation (ARCIx) described in the Methods section. 
Following inspection of the data, we divided the mutants into three 
virulence classes with either a high to moderate likelihood of being 
virulence-attenuated (ARCISlO), a moderate to low likelihood 
(10<ARCI<20) or no likelihood (ARCI>20) (Table S4). The first 
two categories comprised about 40% of the mutants (48 out of 
124) (Table 1). Fifteen of the 48 mutants lacked metabolic genes, 5 
lacked stress response genes, 2 lacked fimbriae/ adhesin genes and 
26 lacked uncharacterized genes. Interestingly, some mutants were 
more outcompeted in the lymph node than in the spleen (or vice- 
versa), suggesting that specific factors are required to colonize each 
organ. To confirm the validity of our hierarchical ranking, the 
individual virulence of a number of mutants from each of the three 
different ARCI classes (11, 5 and 4 mutants, respectively) was 
evaluated in the rat. Eight mutants (variously lacking Pap fimbriae. 
Ail or uncharacterized proteins) were found to be significandy less 
virulent (Table S4) and all had an ARCK20. Furthermore, 65% 
and 20% of the attenuated mutants were from the "high ARCI" 
and "low ARCI" categories, respectively; this suggests the 
existence of a correlation between the ARCI category and loss 
of virulence. 

The attenuation of the mutants described above might 
conceivably result from a secondary mutation generated during 
the construction process. To test this hypothesis, we compared the 
virulence of the wild-type strain and several mutants {Aj>po0988, 
A)ipo2062, \ypo0656, Aypo2560-61, Aypo3369 and Ajpo3991) 
harboring (or not) a wild-type copy of the gene of interest on a 



plasmid (Figure S5). In accordance with the results of the above- 
described infection experiments, all mutants were significandy 
affected in virulence (p<0.05 in a Gehan-Breslow-WUcoxon test). 
The virulence of the complemented mutants Ajipo0988, Aypo2062 
and Ajipo0656 was fully restored. In contrast, the virulence of the 
complemented mutants Aypo2560-61, Aypo3369 and Aypo3391 was 
not restored - indicating that these mutants carry a secondary 
mutation or that the methodology used for complementation was 
inadequate. 

y. pestis requires a limited number of metal acquisition 
systems to counter metal starvation in vivo 

Iron, manganese and zinc play a crucial role in many biological 
processes - including bacterial virulence. Hence, it is not surprising 
that a marked number of metal acquisition systems are up- 
regulated by T. pestis in vivo [7,16]. Notably, Y. pestis overexpresses 
1 2 of the 1 5 proven or potential iron transport systems (Table S5 
and Text SI). It has previously been shown that the absence of the 
Yersiniabactin (Ybt) siderophore-dependent system or the Yfe 
ATP-binding cassette (ABC) iron transporter (but not that of the 
heme transporter Hmu, the hemin storage system Hms or the 
ABC iron transporters Yiu and Yfu) attenuates the virulence of T. 
pestis in a mouse model of bubonic plague [26,27,28,29,30]. Our 
present data (obtained in a different T. pestis strain and a different 
rodent species and strain) confirmed that Hmu, Hms, Yiu and Yfu 
are not required for bubonic plague (Table S5). Furthermore, we 
found that three siderophore-based systems, two ABC ion 
transporters and two iron permeases were not necessary for 
disease production, despite upregulation of the corresponding 
genes in vivo; this is presumably because some of these transport 
systems are not functional in Y. pestis due to frameshift mutations 
or gene disruption by an insertion sequence [31]. 
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Figure 2. Carbohydrate uptake and metabolism by K. pestis in its mammalian host. Genes encoding transporters for fructose (FruAB), 
mannose (ManZYX), maltose (MalEMGKF and MalE2FG), glucose (PtsG), gluconate (GntM) and fructuronate (GntP) and the carbohydrate metabolism 
genes shown in red or purple were all overexpressed in vivo [7]. In virulence tests, mutants lacking a gene highlighted in red were virulence- 
attenuated, whereas those lacking a gene highlighted in purple were not. Genes in grey were not upregulated in vivo and their role in virulence has 
yet to be elucidated; hence, the importance of the metabolic fluxes involving the corresponding enzymes in virulence is unknown (grey arrows). Red 
and purple arrows respectively indicate the metabolic fluxes considered to be important and non-essential for virulence. Yellow arrowheads 
represent the metabolic pathways that glucose and gluconate might follow. 
doi:1 0.1 371 /journal.ppat.1 004029.g002 



In vivo, T. pestis also overexpresses the ZnuABC zinc transporter 
and the MntH manganese transporter (which are also the 
predominant zinc and manganese importers in vitro). However, 
mutants lacking Znu and MntH were found to be fully virulent 
(Figure 1 and Table S5), which agrees with recent findings in a 
different T. pestis strain and a different animal model [32,33]. 

Yersinia pestis requires carbohydrate metabolism when 
colonizing the mammalian host 

Inspection of the expression patterns for genes involved in 
metabolic pathways has suggested that carbohydrates constitute T. 
pestis' main carbon and energy source during successful colonization 
of the mammalian host [7] . The results of our virulence tests indicated 
that the uptake of glucose, gluconate and, to a lesser extent, maltose 
(but not that of fructose, mannose or fructuronate) was important for 
plague production (Figure 2 and Table S3 and S4). Glucose is known 
to be metabolized via glycolysis; however, the latter did not appear to 
be essential for plague production because deletion of the first two 
enzymes in this pathway (Pgi and PfkA) did not affect T. pestis' 
virulence. Although the presence of a complete glycolysis pathway 
was not essential for host colonization, its terminal part appeared to 
be essential because a gpniA mutant was completely outcompeted in 
vivo. The end product of glycolysis (pyruvate) is an important source of 
acetyl CoA, some of which enters the tricarboxylic acid (TCA) cycle 
via reaction with citrate synthase. However, this step (and perhaps the 
TCA cycle as a whole) was not essential for virulence. Thus, T. pedis 
relies on carbohydrate metabolism to produce disease but 



carbohydrates such as glucose are not necessarily directed towards 
the expected metabolic pathways. 

A limited number of genes enable Y. pestis to resist 
nitrosative stress in the bubo 

During infection, Y. pestis is in close contact with PMNs that 
produce inducible nitric oxide synthase. The bacterium uses the 
flavoglobin Hmp to detoxify reactive nitrogen species (RNS) 
released by the PMNs [7]. Of the T. pestis genes that are up- 
regulated in the rat bubo and for which homologs in other 
bacterial species are involved in the response to RNS, aceE that we 
found to be necessary for virulence was also found to be required 
for optimal bacterial growth in vitro (Table S6 and Text SI). 
Therefore, one can hypothesize that aceE is essential for virulence 
but is not necessarily required to counter the host PMNs' 
production of nitric oxide. Of the 8 other loci thought to protect 
T. pestis against RNS in vivo, only nrdHIEF (the ribonucleotide 
reductase operon) appeared to contribute moderately to virulence; 
deletion of this operon was associated with a slightly longer time to 
disease onset (Table S6). 

Further evidence to suggest that Y. pestis may not 
encounter oxidative stress in vivo or may require a 
limited number of genes to resist oxidative stress 

The absence of induction in the rat bubo of (i) the OxyR 
oxidative stress regulator regulon and (ii) the genes known to 
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Figure 3. Incidence of plague in immunocompetent rats (A), immunocompetent mice (B), and neutropenic mice (C) injected 
intradermally with ~20 WT Y. pestis {white circles), \ypmt7. 66c {b\ack circles) or complemented mutant (grey circles). The mutant was 
significantly (p<0.05) less virulent than the WT strain in immunocompetent rats (A), immunocompetent mice (B and C) and neutropenic mice (C). The 
survival curves of mice infected w\th the complemented mutant harboring the ypmtl.66c gene on a high copy number plasmid (pCR2) and the wild- 
type (B) were significantly different (p<0.05). The virulence of Y. pestis lacking ympt1.66c was no greater in neutropenic mice than it was in 
immunocompetent mice (dashed lines) (p>0.05). Data were obtained from groups of 8 or 9 animals. P values were determined using the Gehan- 
Breslow-Wilcoxon test. 
doi:1 0.1 371 /journal.ppat.1 004029.g003 



protect bacteria against peroxynitrite (i.e. the product of the 
reaction between hydrogen peroxide and nitric oxide) suggested 
that T. pestis is not exposed to reactive oxygen species (ROS) in the 
mammalian host [7]. However, a limited number of genes 
involved in bacterial resistance against ROS are known to be 
activated in the rat bubo, and they may encode the major 
detoxifying proteins. Furthermore, several genes involved in the 
detoxification of ROS (which are downregulated or at least not 
upregulated in the rat bubo) were shown to be upregulated in T. 
pestis replicating inside a macrophage-cell line in vitro [34]; 
macrophage infection is thought to be important in the initial 
stages of infection [6]. Of the genes known to protect bacteria 
against oxidative stress, T. pestis overexpresses mntH, fumC, yfiD 
IgrcA) and ibpA in vivo. As mentioned above, MntH is required for 
manganese import (which in turn is involved in protection against 
ROS), whereas FumC and GrcA rescue enzyme reactions when 
FumA and PflB are damaged by ROS [35,36,37]. IbpA protects 
several metabolic enzymes against oxidation [38] . Deletion of the 
ibpA gene was associated with attenuation of virulence because the 
mutant was outcompeted in per-pool infection. However, IbpA is 
also involved in resistance to other forms of stress, such as high 
temperatures [38,39] . In contrast, deletion of fumC, yfiD or mntH 
did not affect the virulence of T. pestis (Table S6). Thus, our data 
provide further evidence for the hypothesis whereby ROS 
production may not be a significant host defense factor in plague. 

Horizontal genetic transfer had an important role in 
acquisition of V. pestis capacity for intracellular survival 

Some of the Y. pestis genes overexpressed in the rat bubo and 
required for virulence in the present study have been acquired 
horizontally. These genes include the P4-like prophage genes 
{ypol089 to ypol 091) and the yprntl. 66c gene (which encodes a 
putative DNA binding protein and is located on the pMTl 



plasmid). A A)ipmtl.66c mutant was outcompeted in a per-pool 
assay (Table S4). Furthermore, the absence of the jipmtl. 66c gene 
decreased the incidence of plague in both rats and mice (Figure 3). 
The virulence of the mutant expressing a WT copy ofj>pmtl.66c 
under the control of its own promoter from a high copy number 
plasmid (pCR2) was almost identical to that of the wild-type strain. 
The complemented mutant killed 100% of animals but caused 
fatal plague with a slight lag (a median survival time of 4 days, vs. 
3.5 for the WT). Even though this slight difference in virulence was 
statistically significant (P = 0.038), this result indicated that loss of 
virulence was certainly due to deletion oiypmtl.66c. Interestingly, 
animals which survived inoculation with the Kypmtl.66c mutant 
had enlarged, purulent lymph nodes but not severe splenomegaly - 
suggesting that the mutant was impaired in its ability to initiate 
skin colonization and/ or to escape from the draining lymph node. 
Consistently, the times to colonization of the lymph node, blood 
and spleen were significandy longer in rats infected with the 
mutant than in rats infected with the WT strain (Figure 4). Taken 
as a whole, our data indicate that a mutant's virulence is 
attenuated because the host manages to contain the infection 
during the initial skin colonization step, for which intracellular 
replication within macrophages is thought to be important [6] . In 
accordance with this model, the A)i/)mil.66c mutant was more 
susceptible to macrophages than the WT and the complemented 
mutant (Figxire 5A). The mutant's failure to initiate infection at the 
fleabite site could also result from its inability to use available 
nutrients or to survive contact with PMNs once released from 
lysing macrophages. The mutant showed a slight but statistically 
significant growth defect in serum (used here as model of nutrient 
acquisition in the skin), when compared with the WT and the 
complemented strain (Figure 5B). Furthermore, the mutant's 
respective levels of virulence in neutropenic and immunocompe- 
tent mice did not differ significandy (Figure 3C); this indicates that 



PLOS Pathogens | www.plospathogens.org 



8 



IVIarch2014 | Volume 10 | Issue 3 | el 004029 



Adaptive Responses of Y. pestis during Plague 



Lymph node 



o 



9-1 

7- 



= 5 



3- 



ND 
I C03- 



Blood 



O 0> 

o a 5- 



Spleen 



hours after challenge 



Figure 4. The change over time in colonization and bacterial load in rats inoculated intradermally with ~20 WT V. pestis (white 
circles) or A/pmfl .66c (black circles). N.D: not determined. Each circle represents one animal. The times to colonization were significantly longer 
in rats infected with the mutant than in rats infected with the WT strain (p<0.05). P values were determined using 2-way analysis of variance. 
doi:1 0.1 371 /journal.ppat.1004029.g004 



YPMT 1.66c is involved in a biological process situated upstream 
of interaction with PMNs. Taken as a whole, our data suggest that 
YPMT1.66c's role in survival within macrophages is important in 
the production of plague. 

Stress, metabolic and uncharacterized genes are required 
for resistance and/or growth in serum 

Approximately 40% of the mutants found to be virulence- 
attenuated in individual tests and classified in the high-to-low 
likelihood class for virulence attenuation (i.e. ARCISIO and 10< 
ARCK20) are known to encode proteins related to metabolic 
pathways. The remaining 60% of mutants lack factors considered 
to be involved in stress resistance (e.g. IbpA or RseC), bacterial 
attachment or serum resistance (e.g. Ail). Indeed, most of the 
remaining mutants are involved in hypothetical or uncharacter- 
ized processes. To determine whether the selected mutants were 
impaired for bacterial growth in vivo rather than for resistance to 
the host immune system, we measured their ability to grow in fresh 
serum (Figure 6). However, the mutants with a 10<ARCI<20 
were not evaluated because we expected few strains in this ARCI 
category to be truly virulence-affected (see discussion) and 
therefore expected even fewer mutants to be serum-sensitive. 
Furthermore, the AaceEF mutant (identified in an individual 
virulence test) was not evaluated because its serious growth defect 



in LB led us to expect a growth defect in serum. Lastly, the 
ARCI^IO mutants lacking genes that encode fimbriae or genes 
found to be non-essential for plague in individual tests were not 
evaluated because (i) we did not expect fimbriae to be involved in 
serum resistance and (ii) mutants that were found to be attenuated 
in per-pool screening but not in individual tests were considered to 
be false positives. Hence, a total of 23 mutants were evaluated for 
serum resistance. Ten of them differed significandy from the wild- 
type strain. In particular, the ArseC, Ajpo0337 and AgpmA mutants 
were sensitive to serum (presumably because of the bactericidal 
effect of the complement), and the AibpA, Aypo3369, Aypo0988, 
Aamn, Aypo0617-0618, Aypo2586-2587 and Aypo0426 mutaiAi had 
significant a growth defect (either because they confer a metabolic 
advantage or effective complement tolerance). 

Discussion 

Here, we presented a novel experimental approach for 
establishing a hierarchical table of Y. pestis genes previously 
identified (in a microarray analysis) as being up-regulated in vivo. 
The genes' ranking in terms of their importance in virulence 
should provide new insights into the pathogenesis of plague. Our 
approach consisted of (i) testing the virulence of individual mutants 
each lacking a locus likely to encode factors important for plague 
pathogenesis and (ii) per-pool screening of a library of mutants 
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Figure 5. YPMT1 .66c is required for intracellular survival in macrophages (A) and for optimal growth in serum (B). Shown are the data 
obtained with the wild-type (white squares and bars), the Aypmtl.66c (grey squares and bars) and the complemented mutant (black squares and 
bars) strains. Data are quoted as the mean (SD) from three independent experiments using macrophages (A) and 5 independent experiments using 
serum from 5 different healthy donors (B). The mutant's survival (A) and growth curves (B) differed significantly from those of the WT and the 
complemented strain (p<0.05 in a two-way analysis of variance). The complemented mutant's and the WT's survival (A) and growth curves (B) did not 
differ significantly (p>0.05 in a two-way analysis of variance). 
doi:1 0.1 371 /journal.ppat.1 004029.g005 
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Figure 6. Stress, metabolic and uncKiaracterized genes are needed for resistance and/or growth in serum. The ability of selected 
mutants to grow in fresh serum was evaluated after 6 and 21 hours. Data are representative of three independent experiments using sera from three 
unrelated, healthy donors. Only the growth curves of the ArseC, Aypo0337, AgpmA, AibpA, Aypo3369, Aypo0988, Aamn, Aypo0617-0618, Aypo2586- 
2587 and Aypo0426 mutants differed significantly from that of the wild-type strain (p<0.05 in a two-way analysis of variance). 
doi:10.1371/journal.ppat.l004029.g006 



composed of mutants strains each lacking a locus with no obvious 
role in pathogenesis and of which ~60% of the strains lacked 
genes that were not directly involved in any known metabolic 
pathway. Our per-pool screening metliod is based on signature- 
tagged mutagenesis (STM) [40] and is closely related to "array- 
based analysis of cistron under selection" (ABACUS) technology 
[41]. We first generated a library of non-polar, tagged deletion 
mutants by using lambda-Red technology [13] (as in ABACUS) 
rather than a transposon mutant library (as in STM). Each mutant 
in the pool was tagged with a unique antibiotic resistance cassette, 
which ultimately allowed the numbers of viable, tagged mutants in 
infected tissues to be quantified by CFU counting on selective 
media. This aspect contrasts with the ABACUS method, in which 
(i) a T7 RNA polymerase promoter is inserted into the target gene 
and (ii) mutants are monitored with a microarray using labeled 
RNAs synthesized in vitro from the T7 promoter and mutant DNA 
that has been fragmented, poly(A)-tailed and amplified prior to 
transcription. Hence, in contrast to ABACUS and other method- 
ologies based on microarrays or deep sequencing (such as STM- 
like or transposon-directed insertion-site sequencing), our method 
does not require nucleic acid purification, RNA synthesis, PCRs, 
microarray hybridization or sequencing [41,42,43]. It therefore 



minimizes the likelihood of introducing bias related to sample 
preparation, hybridization and/or sequencing. One could argue 
that our "low complexity", per-pool approach is limited by the 
number of antibiotic tags available. However, various fluorescent 
tags (such as green or red fluorescent protein), bioluminescent tags 
(blue or red) and enzymatic tags (such as LacZ) can be used to 
increase the number of mutants in the pool. Hence, it would be 
easy to match the pool to the desired complexity and, ultimately, 
to the most relevant infection conditions. In fact, the tag approach 
enabled us to generate a low-complexity pool that corresponded to 
inoculation of a physiological dose of bacteria, i.e. that delivered 
by a flea during a blood meal. 

By injecting pools of mutants under the most relevant infection 
conditions, we confirmed previous suggestions whereby T. pestis 
encounters a bottleneck during dissemination from the skin (Figure 
S4) [20,44,45]. This dissemination bottleneck has been observed 
for several pathogens (including Bacillus anthracis and Fmncisella 
tularensis) and was previously considered to prevent the use of per- 
pool screening methodologies other than when very high doses 
were inoculated [25,41]. Our results do not support the latter 
hypothesis. Our calculation of the ARCI enabled us to address the 
bottleneck problem and thus draw up a hierarchical table of T. 
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pestis genes in terms of their importance in plague. We then 
validated this hierarchical table by testing the virulence of several 
mutants individually. With the exclusion of ail (which was 
identified as needed for plague pathogenesis by other researchers 
while we were performing our experiments [46]), 16 new loci were 
found here to be important for \'irulf;ncc in mice and/ or rats 
(Tables S3 and S4). Among them, the uncharacterizedj)'^mii.6'& 
and ypo0988 genes were required for optimal growth in serum 
(Figures 5B and 6). Furthermore, the uncharacteriz(-d ypo()337 
gene (previously identified by signature-tagged mutagenc'sis in 
Tersinia pseudotuhenulosis [45], the recent ancestor of 7. pestis [47]) 
was essential for serum resistance (Figure 6); this finding suggested 
that YPO0337 is a new, uncharacterized virulence factor. Thirty- 
two mutants identified as being virulence-attenuated in per-pool 
screening were not reassessed individually in the present study 
(Table S4). Eleven and 20 of these mutants belonged to the 
ARGILS 10 and ARCK20 classes respectively. Based on the 
virulence of individual mutants in the different ARCI classes, we 
estimate that 65% of the mutants with an ARCKIO and 20% of 
those with an ARCK20 might truly be virulence-attenuated. 
Hence, at least 12 other loci can conceivably be added to the 
current list of 16 new loci found to be important in plague 
pathogenesis. Accordingly, several genes were required for survival 
or optimal growth in serum (Figure 6). However, we cannot rule 
out the possibility that some of the mutants identified here are 
defective in virulence because they harbor a secondary mutation 
(Figure S5). The apparent absence of complementation might also 
result from use of an inadequate method or other, uncharacterized 
features. 

Our virulence testing results suggest that in the mammalian 
host, T. pestis (i) relies on carbohydrates (i.e. glucose, galactans, 
gluconate and perhaps maltose) as its carbon source and (ii) shifts 
to anaerobic respiration (maybe using dimethyl sulfoxide (DMSO) 
or glycerol) and/ or fermentation (Figure 2 and Tables S4, S5, S6). 
The data indicated that glucose (imported via the PtsG 
transporter) and galactans were channeled towards the synthesis 
of UDP-glucose rather than to glycolysis. Firstly, Y. pestis does not 
encode a functional glucose-6-P dehydrogenase (Zwf, the enzyme 
that catalyzes the first step in the metabolism of hexose via the 
pentose phosphate pathway [48]. Secondly, deletion of the first 
two enzymes in the glycolysis pathway (Pgi and PlkA) did not 
impact virulence (Figure 2 and Table S3). Whereas carbohydrates 
do not feed into the upper part of glycolysis, gluconate (taken by 
the GntM transporter) can potentially supply glyceraldehyde-3- 
phosphate (G3P) to the lower part of the pathway. The latter 
metabolite appeared to be essential for plague production, since a 
gprriA mutant was completely outcompeted in vivo (Figure 2 and 
Table S4). Glyceraldehyde-3-phosphate is an important carbon 
source in hexose and amino-acid anabolism and pyruvate 
production [49]. However, we do not ha\e any e\idence to 
suggest that G3P is involved in hexose and amino-acids anabolism 
in Y. pestis. In particular, SerA (YPO0914, which catalyzes the first 
step in the biosynthesis of serine from glycerate-3P that is 
produced from G3P) was not necessary for plague. Whereas 
G3P's role in hexose and amino-acid anabolism in Y. pestis has yet 
to be confirmed, this compound is likely to be crucial for the 
production of pyruvate in vivo because a mutant lacking the 
pyruvate dehydrogenase AceEF was unable to produce fatal 
plague when injected at a dose known to kill most animals. Lastly, 
the facts that (i) the tricarboxylic acid cycle (TCA) is down- 
regulated in vivo [7], (ii) the mutant lacking the acotinase A (part of 
the TCA cycle) is fully virulent (unpublished data) and (iii) the 
mutant lacking the heat-stable iron-independent fumarate enzyme 
(FumG, part of the TCA cycle) was not outcompeted in per-pool 



virulence screening in contrast to the mutants lacking DMSO 
reductase (DmsABC, YP02965-2969) and glycerol-3P dehydro- 
genase (GlpABC) - both involved in anaerobic respiration - 
suggested that Y. pestis relies on anaerobic respiration to produce 

plague. 

Our virulence testing results also suggest that in the mammalian 
host, Y. pestis experiences iron, RNS and temperature stresses but 
not ROS stress (Figure 2 and Tables S4, S5, S6). This supports the 
hypothesizes generated on the basis of in vivo gene expression 
profiling [7] . However, we tested a small number of genes involved 
in bacterial resistance against ROS (Table S6) because few 
oxidative stress response genes are activated in the rat bubo [7]. 
Several of the genes involved in the detoxification of ROS which 
are not upregulated or even downregulated in vivo are upregulated 
in Y. pestis replicating inside macrophage-cell fine [34] - a step 
suggested to be important in the initial stages of infection [6]. The 
difference in expression levels found in rat bubo and in 
matTophages cell-line in vitro may indicate that Y. pestis encounters 
ROS during the initial stage of infection but not later on. 
However, the detoxifying hydroperoxidase katG gene (the expres- 
sion of which is upregulated ~10 fold during replication within the 
macrophage) was found to be dispensable for plague pathogenesis 
[50]. Hence, the role of ROS during plague infection remains to 
be accurately established by systematic deletion of the ROS 
resistance/responding genes and characterization of the resulting 
mutants. 

Our screening results also suggested that Y. pestis only needs a 
limited number of genes to successfully infect its mammalian host. 
This might be attributable to redundant genes that constitute a 
"spare wheel" system. For instance, serine biosynthesis enzyme 
SerA (YPO0914) was not required for virulence but its loss could 
be rescued by a second putative SerA (YP01288) enzyme;. 
Alternatively, the requirement for a limited number of genes may 
be because only one of several genes with overlapping functions 
are needed in a given context. For example, Y. pestis upregulates 
both idnK &nd encoding gluconate kinases but only gntV 
required for virulence (Figure 2) [51]. One can hypothesize that 
IndK is needed to establish an infection in the flea. Lastly, the 
requirement for just a limited number of genes might indicate 
ancestral needs for colonizing other environments. Hence, one 
would expect at least some of the redundant genes to be 
inactivated and/ or deleted from the genome as the bacterium 
retains only the most appropriate system for its current lifestyle. 
This is consistent with the suggestion that Y. pestis is undergoing 
genome decay [52]. Data on the role of iron acquisition systems in 
plague provide support for this decay scenario. Of the four 
siderophore-based uptake systems overexpressed by Y. pestis in vivo 
(yersiniabactin, pseudochelin, aerobactin and alcaligin) and whose 
homologs in other bacteria are known to have a role in lifestyle, 
only the Yersiniabactin system is crucial for plague - presumably 
because the aerobactin and alcaligin systems are not functional in 
Y. pestis (Table S5) [7,26,53,54]. 

In addition to genome decay, the emergence of Y. pestis has been 
characterized by the horizontal acquisition of the plasmids pPst 
andpMTl/pFra [52]. The two plasmids have had a crucial role in 
the emergence of flea-borne plague [55,56,57,58,59,60,61]. pPst 
encodes a plasminogen activator that is thought to prevent the host 
from containing the Y. pestis infection at the dermal fleabite site. 
pFra encodes the pseudocapsule Fl, whose antiphagocytic activity 
and/or ability to interfere with interleukin- 1 signaling increases 
the incidence of plague. pFra harbors also several uncharacterized 
regions (located near the caf operon and elsewhere) which may 
encode virulence factors [62] . None of these regions contained the 
uncharacterized gene ypmtl.66c (encoding a protein that is 
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conserved among T. pestis strains). Here, we found that that 
ypmt 1.66c is also important for maximum disease production at the 
dermal infection site (Figures 3 and 4). Our observation provides a 
new example of the importance of horizontal gene acquisition in 
the emergence of flea-borne plague. 

The deletion oi ypmtl.6(k reducc-d Y. pestis' ability to resist 
macrophages, and sUghtiy hampers the bacterium's ability to grow 
in normal serum (Figure 5) - suggesting thai ymptl .66c is partially 
needed for resisting complement and/or to acquiring nutrients 
from the serum (and perhaps also at the fleabite site). We did not 
test whether YPMT 1.66c is needed to confront the PMNs that are 
recruited en masse in infected tissues. Even if this is the ca.se, 
YPMT1.66c's role as an intracellular survival factor appears to be 
a prerequisite for resistance to PMNs because the virulence of a 
mutant lacking YPMTl.SGc was not restored in neutropenic mice 
(Figure 3C). Taken as a whole, our data provide support for the 
scenario in which Y. pestis regurgitated at the fleabite site is taken 
up by macrophages and indicate that YPMT 166c is required for 
replication within the phagocyte. This scenario is consistent with 
the current model of plague, in which the conditions for 
intracellular replication at the initial stage of infection (i.e. at the 
dermal flea bite site) are critical for the subsequent production of a 
bubonic infection [6,9]. 

Conclusion 

We leveraged our previously collected in vivo gene expression 
profiling data on Y. pestis [7] to generate and screen a non-polar 
tagged deletion mutant library composed of bacteria lacking genes 
known to be upregulated in vim. This approach constituted the first 
integrated attempt to identify the molecular mechanisms of the 
pathogenesis of plague in the most relevant biological model of 
infection - the mammalian host. Our approach could be 
transposed to other pathogens and notably those for which a 
dissemination botdeneck has been described [24,25]. Our data 
provided further insights into the pathogenesis of plague by 
identifying new genes important for plague pathogenesis (including 
several previously uncharacterized sequences). The fact that these 
genes are shared with several other pathogens suggests that they 
belong to a core set of bacterial pathogenesis genes. Confirmation 
of these virulence results (by individually testing the virulence of 
mutants selected from the pools and by performing complemen- 
tation experiments) and further characterization of these new 
identified genes wiU be necessary. Future studies should therefore 
improve our understanding of pathogenesis and, in particular, 
yield putative antibiotic targets for several diseases. 

Supporting Information 

Figure SI Equivalent incidences of plague in Brown-Norway 
rats inoculated intradermaUy with either ~IQ Y. pestis C092 or 
-10 Y. pestis C092 mini-Tn7T;.-.^eol (/)>.045 in a Gehan- 
Breslow-Wilcoxon test). Data were obtained from one experiment 
on a group of 8 animals. 
(TIF) 

Figure S2 Incidence of plague in Brown-Norway rats inoculated 
intradermaUy with an input pool of five different Y. pestis mutants 
(20 CPUs per mutant and 100 CFUs in total) (A-F). Groups of 10 
animals were infected with each pool and four groups were 
infected on the same day. 
(TIF) 

Figure S3 Bacterial loads in (A) the lymph nodes draining the 
inoculation site and (B) the spleens from the rats used in the 



per-pool screening. Rats were inoculated intradermaUy with an 
inoculum comprising five mutants (20 CFUs per mutant and 
100 CFUs in total). The bacterial loads recovered from the lymph 
nodes and the spleens from the various groups of rats did not differ 
significantiy (p<0.05 in a one-way analysis of variance with 
Bonferroni correction for multiple comparisons). 
(TIF) 

Figure S4 A bottleneck affects the dissemination of Y. pestis from 
the intradermal inoculation site. The figure shows the bacterial 
loads recovered from the lymph nodes and the spleens from a 
group of 10 rats (Rl to RIO) inoculated intradermaUy with five 
mutants (20 CFUs per mutant and 100 CFUs in total). TypicaUy, 
each animal was predominanfly colonized by just one of the pool's 
mutants. Red s)rmbols indicate the predominant mutant strain. 
Grey symbols indicate mutants with at least 30% of the CFUs of 
the predominant strain. 
(TIF) 

Figure S5 Incidence of plague in mice infected with Y. pestis 
lacking selected uncharacterized genes. (A-F) Animals were 
injected intradermaUy with ~10 CFU of Y. pestis WT (white 
circles), the mutant lacking the gene of interest (white squares or 
white triangles) or the mutant harboring a wild-t}pe copy of the 
gene of interest (grey squares) under the control of a pAra 
promoter (A) or its own putative promoter (B-F). AU six mutants 
were significantiy affected in virulence compared to the wild-type 
strain (p<0.05). In contrast to the complemented hypo3369, 
Aypo2560 -2561 and A)>po3991 mutants, the virulence levels of the 
Aypo0988, Aypo2062 and Aypo0656 mutants harboring a wUd-type 
copy of the gene of interest ("c") differed significantly (p<0.05) 
from that of the original mutant strain but were similar to that of 
the wUd-type strain. . Groups of eight animals were used. 
(TIF) 
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